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Abstract 



Dark energy of phantom or quintessence nature with an equation of state parameter w 

almost equal to — 1 often leads the universe evolution to a finite-time future singularity. An 

o ■ elegant solution to this problem has been recently proposed [Ij under the form of the so- 

O I called Little Rip cosmology which appears to be a realistic alternative to the ACDM model. 

A viscous Little Rip cosmology is here proposed. Whereas generically bulk viscosity tends 

to promote the Big Rip, we find that there are a number of situations where this is not the 

case and where the formalism nicely adjusts itself to the Little Rip scenario. We prove, in 

CN| I particular, that a viscous fluid (or, equivalently, one with an inhomogeneous (imperfect) 

''^ ' equation of state) is perfectly able to produce a Little Rip cosmology as a purely viscosity 

.^ . effect. The possibility of its induction as a combined result of viscosity and a general 

I I I (power-like) equation of state is also investigated in detail. To finish, a physical, inertial 

Q ■ force interpretation of the dissolution of bound structures in the Little Rip cosmology is 

presented. 

1. Introduction. — Viscous cosmology models have been increasingly popular lately. From a 
^ . hydrodynamical viewpoint this is quite a natural development, as the inclusion of the viscosity 

coefficients (there are two of them, shear and bulk) means physically that one departs from the 
case of an ideal fluid and incorporates the deviations from thermal equilibrium to flrst order. 
The case of an ideal (nonviscous) fluid is, after all, quite an idealized model, useful in practice 
for many situations but not so for all, especially not when fluid motion near solid boundaries 
is involved. Also under boundary-free conditions-such as in free turbulence — viscosities are 
physically most important. In a cosmological context, as the cosmic fluid is assumed to be 
spatially isotropic the shear viscosity is usually left out. Anisotropic deviations such as those 
encountered in the Kasner universe are quickly smoothed out being the only coefficient left (, 
the bulk viscosity. Some of the early treatises on viscous cosmology are [2]^[5]. Recently, work 
has aimed at the study of the dark energy universe when the dark fluid is a viscous one [6]. 
In [7] a dark fluid with a time dependent bulk viscosity was considered as a fluid obeying an 
inhomogeneous (imperfect) equation of state. 
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One place where the appearance of bulk viscosity in the cosmic fluid should be expected to 
play an important role, is in the Big Rip (BR) phenomenon, i.e. the singularity of the universe 
in the future (for pioneering papers on the BR, see [SI IS])- The Big Rip (or type II, III and IV 
Rip |TOl [11]) means that one or more of the physical quantities go to infinity at a finite time t 
in the future. Mathematically this implies divergent integrals, typically met when one uses the 
Friedmann equations to express t as an integral over the density p. For phantom/quintessence 
dark energy models which develop a finite-time future singularity of one of the four known 
types pjj, the following scenarios to avoid future singularities were proposed: quantum effects 
[12], modified gravity [13], and coupling with dark matter [H]. 

Recently, a novel scenario has been proposed in [Ij, the so-called "Little Rip" (LR). Models 
are there examined in which the nonviscous dark energy density increases with time (equation 
of state parameter w < —1), but w — )■ — 1 asymptotically, in that way avoiding the future 
singularity. Typically, this was found to occur when the scale factor increases rapidly with 
time, like a{t) ~ exp(exp(t)) or higher exponentials. And this brings us to the topic of the 
present paper, namely to examine the consequences of endowing the fluid with a bulk viscosity. 
In most cases, the presence of a bulk viscosity tends to promote the Big Rip phenomenon; an 
example being presented at the end of section 2. However, it turns out that the formalism can 
nicely fit the Little Rip scenario under special conditions. One such example will be considered 
in detail in section 4. Particular underlying conditions for this case to occur will be thoroughly 
discussed for a large class of models. Specific examples are when the bulk viscosity is constant, or 
is inversely proportional to the Hubble parameter. A mixed effect of viscosity and a generalized 
(power-like) equation of state in the non-viscous sector can also lead to a LR cosmology. It 
will be explicitly demonstrated that, for the standard equation of state p = wp, just taking 
into account the special viscous sector may bring the universe evolution to a LR. To finish, an 
inertial force interpretation of the dissolution of bound structures in the LR cosmology will be 
developed. 

We may here notice that there is in general a problem with a possible transient acceleration 
in the universe: the acceleration may turn back to a usual deceleration regime in the future. 
See, for instance, the papers [TJ] discussing this theme. In principle, one may expect such 
scenarios also in modified gravity. In our present model, however, the Little Rip cosmology 
turns out to be stable (effectively the phantom-like regime), so we do not expect such a transient 
acceleration behavior in the model under consideration. 

2. Basic formalism: General viscosity. — We assume first that that the cosmic fluid possesses a 
shear viscosity r] in addition to the bulk viscosity (. We put ks = c = 1, take the Minkowski 
metric in the form (-+++), sum Greek indices from to 3, and Latin indices from 1 to 3. 
If U^ = {U^, W) denotes the four-viscosity of the fluid, we have in comoving coordinates 

Let Qf^u be a general metric, and introduce the projection tensor h^^, = g^^, + f/^t/,^.Then the 
rotation and expansion tensors can be expressed, respectively, as 



1 



u 



being the scalar expansion d = 9^ = W^-^. In a FRW space 6 = 3H, with H the Hubble 
parameter. This relationship will be made use of in the following. The third tensor to be 
defined is the shear tensor, 

(^^w = d^jy — Hh^y. (2) 



It is traceless, a'^ = 0. The following decomposition of the covariant derivative U^-i, is often 
useful, 

U^,-^ = ujf,y + a^,y + Hh^y - A^Uy, (3) 

where A^ = [/^ = W^U^-a is the four- acceleration of the fluid. 

Consider now the fluid's energy-momentum tensor T^^, assuming for simplicity that the 
temperature T is constant, 

T^u = pU^U, + {p- KH)Ku - 2w^u. (4) 

Here 

Pes=p- SCH (5) 

is the effective pressure, lower than p because of the inequality C ^ which in turn is a 
consequence of thermodynamics. 

The entropy per unit comoving volume is S* = na, where n is the particle number den- 
sity and a the entropy per particle. This corresponds to an entropy current four-vector and 
corresponding four- divergence 



■,fl rp ^/l!^^ T^ rp 



S^ = naW^, S'^.,, = 4a^,a^^ + -^E\ (6) 



Thus, both viscosity coefficients, ri and C, contribute to the heat production in the fluid. 

We now specialize to the case of an isotropic fluid, ?7 = 0, and introduce the equation of 
state 

P=-P-/(P), (7) 

/(p) being an arbitrary function. What is most important in viscous cosmology, is not p but 
instead the effective pressure peff- We generalize the expression ([5]), assuming that C(-f^) is a 
general function of H. Then, with 'iC^H — > ^(-f/^) we express the effective pressure as 

Peff = -p-/(p)-e(i^), (8) 

which is also called cosmological inhomogeneous (imperfect) equation of state [7] (in general 
one should allow the function ^ to depend on derivatives of B. too). 
Consider now the spatially flat FRW line element 

ds^ = -dt^ ^G?{t)(b^. (9) 

The first and second Friedmann equations are, respectively [H = a/ a), 

H' = y Gp, ^ + Ih' = ^^G[p + /(p) + an)] (10) 

(note that it is essentially the effective pressure that appears on the right hand side). From the 
conservation of energy, T^" -^y = 0, we get 

p-3f{p)H = 3aH)H. (11) 

The special case where ({p) is set proportional to H is of interest: 



((p) = 3tH = T,/24^, (12) 



with T the proportionahty constant. The ansatz (fT2!) should be quite reasonable physically, as 
it implies that the viscosity, i.e. the deviation from thermal equilibrium, is most pronounced 
when the motion of the universe is most violent. In this case, we obtain 

v/24^ypoP'/M/(p)/P + 24vrGr]- ^ ' 

Here the integration starts from to = 0, the present time. The Big Rip is inevitable if this 
integral is finite for p — )■ oo. 

It is instructive to make the ansatz [71 [1] 

/(P) = Ap", (14) 

where A and a are constant. If the fluid is nonviscous, r = 0, it is seen that the value a = 1/2 
represents a borderline between Big Rip and non-Big Rip. If a < 1/2, the singularity is avoided. 
However, for a viscous fluid the situation is different: even if f{p)/p — > for p — )• oo the integral 
reduces effectively to / p~^^'^dp oc p~^/^ for large p. 

We thus arrive at the following important conclusion: ii ( oc H the Big Rip is inevitable, 
regardless of the form of the function / (p) . 

3. Little Rip cosmology: no viscosity. — We first reproduce some essentials from Ref. [Ij. When 
^{H) = it follows that the equation H = p/(3/(p)) is satisfied by the following integral 
expression for the scale factor: 

We have here put the scale factor at present time, Oq, equal to unity. The first Friedmann 
equation then yields 

"" (16) 
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We assume hereafter that /(p) is given by the power law f lT^ . and put a = 1/2. Then 
/(p) = ^a/P, and 

p=-p- A^p. (17) 

From Eq. (IT3l) we obtain in the nonviscous case 

t = ^^iln^, (18) 

^flMG A po 

showing that p — )■ oo is not reached until an infinite time t has elapsed. This is the Little Rip 
phenomenon. Note that it is called this way because it also produces disintegration of bound 
structures, as we discuss latter. 

Using Eq. (fT5|) we can also express p as a function of a. 

Again using the first Friedmann equation we get 



" daja \ 3 

, l + ^lna = V^^^'' ^ ^ 



which upon integration permits us to find the scale factor as a function of time, 



a{t) = exp 



3A 



exp 



3A 
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(21) 



As explained in Ref. [T], this LR cosmology is effectively non-singular, because the singularity 
only occurs in the infinite future. 

4. Little Rip cosmology: constant viscosity function CiH)- — Now, we move on to consider the 
viscosity-dependent governing equations. Friedmann's first equation in (fTOj) remains formally 
unchanged, while the second and (ITT]) become altered. It is seen that the formalism becomes 
easily manageable if the following condition on the viscosity holds. 



^{H) = ^0 = const. 



In that case we obtain from Eq. ( ITT]) 



a = exp 



dp 



'Po3[eo + /(p)] 

and further use of the first Friedmann equation then yields 

1 r dp 



(22) 



(23) 



t 
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(24) 



If we assume f{p) = A^/p as before (it implies that the dimension of A in geometric units is 
"2^ we get 

1 1,1^,^0 + ^yp 



cm 



t 



^/6^A ^o + ^Vpo 
This equation can be inverted to give p as a function of t, 



Pit) 
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(25) 



(26) 



Thus p — )■ 00 can be reached, but the universe needs an infinite time to do so. This is precisely 
the Little Rip phenomenon, now met under viscous conditions. Recall again that this conclusion 
rests upon the equation of state being adopted in the form ( IT71) . We see that p(t) increases more 
rapidly in the presence of viscosity. It may be interesting to compare it with the turbulence 
effect to dark energy [T6] . 

The effective parameter determining the infiuence from viscosity is the nondimensional 
quantity e = ^o/(^a/Po)- A weak infiuence from viscosity, which is probably the case most 
expected, implies that e <^ 1. The nondimensional quantity ^o/(^VpW) will rapidly become 
much less than e as time increases. 

5. More general cases. — The above analysis can be extended to much more general cases, 
concerning the dependencies of both the function /(p) and the bulk viscosity, (, in terms of 
the density p (or Hubble constant H, see Eq. flT^ ). The most general situation which will be 
here discussed is: 

f{p) = Ap^-^'/^ C = hp\ (27) 



Note that the particular case studied in the previous section corresponds, in this notation, to 
7 = —1/2, /3 = 0. Proceeding as before, we now obtain for Eq. 



I f dp 
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This expression can be integrated generically in terms of an hypergeometric function, the result 
being 



, = ^.^£Z^+ ^-' 



y/2A^ { 7&(7 - /3) 7^(7 - /?) . 



P 



7-/3 7-/3 /3 



(29) 



Some particular cases deserve special attention. When the bulk viscosity is constant (that 
is, 7 = 0) and also /3 = (as in the previous section), we obtain again the result that t has a 
logarithmic dependence on p 

t= ^ ^ In^, (30) 

^/2A^{A + b) Po 

and, proceeding as above, the Little Rip phenomenon shows up. Notice that this was obtained 
before for the case when the bulk viscosity was inversely proportional to the Hubble constant 
(or, what is the same, to the square root of the density), while here it is the bulk viscosity itself 
which is strictly constant. In fact, it is easy to see that this result is more general: keeping 
/3 = 0, for any 7 < we have a LR, while for 7 > a BR occurs. It is plain that a similar 
situation appears in the opposite case: when 7 = 0, for /3 < we have a LR, while for /3 > a 
BR. 

In the subclass of cases when 7 = /3 > one obtains that the value p = 00 is reached at a 
finite and positive time, namely, 

tBR = ^°'' , (31) 

V24^(A + 6)7 ^ ^ 

so that we are back to the Big Rip situation in all these cases. Things change, however, when 
7 = /3 < 0, in which circumstance no singularity can be ever reached in finite time. 

An asymptotic analysis of the general formula (!29|) . which encompasses the whole class of 
dependencies here discussed (1271) . yields a explicit result for each particular situation one may 
be interested in, within the three possibilities: absence of a Rip singularity (e.g., 7 > /3, with 
7 > 0), eventual formation of a Little Rip singularity (that is, one that is reached asymptotically 
at infinite time), and formation of a Big Rip singularity at a finite time in the future (e.g., 7 > /3, 
with 7 < 0). The interplay between the dependencies of both the function f{p) and the bulk 
viscosity, (, in terms of the density p (or Hubble constant if), in getting a specific cosmology 
in which one or another of these possibilities is realized, is explicitly manifest in the calculation 
carried out in this section. 

6. On the physics of the viscous LR. — In order to better understand the influence of cosmolog- 
ical viscosity in physical terms, we start again from Eq. ([8]) but now, for simplicity, 

p = wp~ ^H) , (32) 



that is, the ordinary EoS just modified by addition of the viscosity term. The conservation of 
the energy density reads 

= p + 3H[{l + w)p-^iH)] . (33) 



Consider now the following Little Rip model: 

H = Hoe^' 



(34) 



with Hq and A positive constants. The parameter A in [Ij and Eq. ( !T4|) corresponds to \/2\ 
in dHl]) and we know it is bounded as 2.74 x 10"^ Gyr"^ < A < 9.67 x 10"^ Gyr"^ by the 
Supernova Cosmology Project observational data [IT]. In [1], the model ( l34l) gives almost 
identical behavior of the distance versus redshift and, therefore, the simple model flM|) turns 
out to be very realistic. 

By using the first FRW equation, we find 



and therefore 



Then, from Eq. ([MD 






n-^^ H-2^2At _ 6A 2 

P - -^^0^ - -^^ ■ 



= ^H^ + 3H 



K^ 






which can be solved with respect to ^{H) as 






(35) 
(36) 
(37) 

(38) 



Conversely, if we start with a perfect fiuid whose equation of state is given by (132 p and 
there is a solution which realizes the LR as (I34p . The remarkable point is that such LR is 
induced purely by the viscosity sector. 
Instead of flMj) . we may now consider 



H = Hoe^ 



(39) 



Here Hq, C and A are positive constants. It has been shown in pQ that this model fl5^ can be 
realistic too. We find p and p as 



^ = ^^-' 



ny 



p 



6CA 



i/^e^^^'e^* 
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(40) 



while ^{H) has the following form: 

cm 



i(i±i^i,^ + !4ffto^. 



Using the first FRW equation, p = -^H"^, the conservation law (133!) can be rewritten as 



(41) 



= —H + 3 



K^ 



3(l + ti;) 



H'-^H) 



(42) 



In the case of (l38l) . we obtain 

= H-\H, (43) 

whose unique solution is flMl) up to the constant shift of the time coordinate t — )■ t + to- On 
the other hand, in the case of fj4Tl) . we get 

= H~\H\n^. (44) 

-no 

For this Eq. (14^ . besides the solution ( 139|) . the de Sitter solution H = Hq appears, which can 
be regarded as a special case corresponding to C — ;■ in (!39|l . The stability of the de Sitter 
solution can be now investigated by performing the perturbation 

H = Ho + 6H. (45) 

Then, Eq. dS]) yields 

= 5H- X6H . (46) 

We find that 6H grows exponentially 6H oc e^* and, therefore, the de Sitter solution H = Hq 
is unstable. This indicates that the universe could actually evolve into the LR solution (15^ . 

When the universe is expanding, the relative acceleration between two points, whose distance 
is /, is given by Id /a [a is the scale factor). A particle with mass m at a certain point will be 
subject to an inertial force 



inert 



mVd/a = ml(H + H^] , (47) 



whenever it is observed from another point. Let us assume the two point force is bounded by 
the bounding force Fq. If Finert is positive and bigger than Fq, the bound on the two points 
is destroyed. This is a "Rip" caused by the accelerating expansion. The "Rip" always occurs 
when H diverges, which corresponds to a "Big Rip" [8]. Even if H is finite, if if > diverges, 
which corresponds to a Type II or "sudden future" singularity [IOl[TT], the strong rip will always 
occur. Even if H is finite in finite future, in the case when H becomes infinite in the infinite 
future, the inertial force fH7|) becomes larger and larger and destroy any bounded object. This 
has been dubbed as "Little Rip" pQ. In this case and for realistic models, the time ^lr from 
present until the destruction of the Earth-Sun system has been estimated. In the case of (IMjl 
it amounts to ^lr ~ 8 Tyrs and in the case of ( 139|) to ^lr ~ 146 Gyrs. 

In summary, we have here considered the role of a viscous (or inhomogeneous (imperfect) 
equation of state) fluid in a Little Rip cosmology. Despite the earlier observations that viscosity 
basically supports the Big Rip singularity, we have demonstrated that it is also able to give rise 
to a non-singular. Little Rip cosmology, which is considered to be a viable alternative to ACDM 
cosmology. In particular, constant bulk viscosity and a viscosity inversely proportional to the 
Hubble rate have been considered. We have shown that in those cases a Little Rip cosmology 
may naturally emerge. It is remarkable that for a standard fluid, p = wp, the only influence of 
viscous effects can naturally drive the universe to a Little Rip type evolution. What is more, 
in some cases the Little Rip cosmology may be even more stable than de Sitter space itself. 
This invites to a very careful investigation of all properties which are specific of the Little Rip 
cosmology, which avoids the well-known problems of the Big Rip universe. 
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